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Abstract. The present work extends the weU-known thermodynamic relation C = 
{5E'^) for the canonical ensemble. We start from the general situation of the ther- 
modynamic equilibrium between a large but finite system of interest and a general- 
ized thermostat, which we define in the course of the paper. The resulting identity 
(SpSE) = 1 + {SE"^) d^S (E) jdE'^ can account for thermodynamic states with a nega- 
tive heat capacity C < 0; at the same time, it represents a thermodynamic fluctuation 
relation that imposes some restrictions on the determination of the microcanonical 
caloric curve /3 (E) — dS (E) /dE. Finally, we comment briefly on the implications of 
the present result for the development of new Monte Carlo methods and an apparent 
analogy with quantum mechanics. 
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1. Introduction 

From the standard perspective of Statistical Mechanics, it is costumary to start with 
the Gibbs canonical ensemble given by: 

iP, N) = — ^ exp I-PHn] , (1) 



Z{(3,N) 

which provides the macroscopic description of a Hamiltonian system Hn in 
thermodynamic equilibrium with a heat bath (a very large heat reservoir) at constant 
temperature T, where (3 = l/ksT. Hereafter, the Boltzmann constant fcs is set to 1. 
In this context, other thermodynamic parameters of the system, like the volume V or 
an external magnetic field H, are also admissible, but we assume along this paper that 
every parameter remains constant. A straightforward consequence of using this kind 
of statistical ensemble is the relation between the heat capacity C = dE/dT and the 
canonical average of the energy fluctuation {6E'^): 

C = f3' {6E') , (2) 



Thermodynamic fluctuation relation for temperature and energy 



2 



which leads to the non-negative character of the heat capacity within the canonical 
description, e.g. Cy > for a fluid with volume V, or Ch > for a magnetic system 
in external magnetic field H. This same result can be also derived from the stability 
analysis in the framework of the standard Thermodynamics, therefore, it is usually 
claimed in several classical textbooks that those macrostates, in which this condition 
is not satisfied, are thermodynamically unstable and cannot exist in Nature (see, for 
instance, in the section §21 of the Landau & Lifshitz book [Ij). 

Surprisingly, macrostates with negative heat capacities have been actually observed 
in several systems belonging to different physical scenarios. For example, in the 
mesoscopic short-range interacting systems like the nuclear, atomic and molecular 
clusters [21 l SI [5], as well as the long-range interacting systems like astrophysical 
ones P, El [8], all of them often referred as nonextensive systems^. This observation 
illustrates the limited validity of certain standard results of classical Thermodynamics 
and Statistical Mechanics in such contexts [9]. We are going to provide two illuminating 
examples below, but first, let us explain how negative heat capacities arise in the 
thermodynamic description. 

The usual way to access macrostates with negative heat capacities is by means of 
the microcanonical description. The fundamental key is to rephrase the heat capacity 
C = dE/dT in terms of the Boltzmann entropy S = \nW. Starting from the definition 
of the microcanonical inverse temperature of the system: 

we obtain then the second derivative of the entropy: 

dE^~ T^dE^ ~ ydEj \dE^) ■ 

Since the Boltzmann entropy is a geometrical measure of the microcanonical ensemble, 
it demands neither the extensive and concave properties usually attributed to its 
probabilistic interpretation: 

Se[p] = -^Pk^npk, (5) 

k 

nor the consideration of the thermodynamic limit. Thus, the expression (jl]) states that 
negative heat capacities are directly related to the presence of macrostates with a convex 
microcanonical entropy, d'^S/dE'^ > 0. 

To our knowledge, the existence of negative heat capacities was firstly pointed out 
by Lynden-Bell in the astrophysical context in seminal papers [SI El- Interestingly, the 
presence of this anomaly plays a fundamental role in understanding the evolution of 
these remarkable physical systems [10]. A simple astrophysical model that shows an 

I Roughly speaking, a system is nonextensive when it cannot be trivially divided in independent 
subsystems, which is explained by the existence of underlying interactions or correlation effects whose 
characteristic length is comparable or larger than the system linear size. Thus, the total energy in such 
systems is nonadditive, and frequently, they are spatially nonhomogeneous. 
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Figure 1. Microcanonical caloric curve of tlic Antonov isothermal model (after [TT]). 
The branch AB represents macrostates where the heat capacity Cy < 0. Notice 
also the nonexistence of equilibrium macrostates when T < Tb — OAGM/R or 
E < Ea — — 0.335GM^/i?, which is related to the occurrence of a gravitational 
collapse at low temperatures (canonical ensemble) or at low energies (microcanonical 
ensemble) . 

energetic region with Cy < is the Antonov isothermal model [11]: a system of self- 
gravitating identical point particles with a total mass M enclosed in a rigid spherical 
container of radius i?, whose microcanonical caloric curve is depicted in FIGiH 

The existence of macrostates with negative heat capacities can be also related 
to the occurrence of phase coexistence phenomenon or first-order transitions in finite 
short-range interacting systems prior to any thermodynamic limit |12l [13], HH [15] . Such 
a relationship is clearly illustrated in FIGE] for the g = 10 states Potts model on a 
square lattice N = L x L: 

Hn = J2{1- 5...,) , (6) 

with periodic boundary conditions, where the sum considers only the nearest neighbor 
interactions [16j. Here, the microcanonical caloric curve OS /BE in terms of the energy 
per particle e = E/N shows a backbending indicating the presence of an anomalous 
region where Ch < at if = |H| = c(§|. In the neighborhood of the critical point 
f3c — 1-42, the bimodal character of the energy distribution function within the canonical 

§ The case with |H| 7^ is described by the Hamiltonian Hn — Hn — H • M, where the total 
magnetization M — and = [cos {na-i) , sin (ncri)] with k =^ 27r/q. 
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Figure 2. The energy distribution functions within the canonical ensemble, obtained 
from Eq.(l7]) by using Vl{E) estimated from Monte Carlo calculations (see sec. 13.11 
below), indicate the occurrence of phase coexistence phenomenon in the g = 10 states 
Potts model on the square lattice N = L x L with periodic boundary conditions 
(e = E/N). Notice that the microcanonical states with negative heat capacities are 
practically inaccessible by using a thermostat with constant temperature, since these 
macrostates are canonically unstable. 

ensemble: 

Pc (E) dE = ^ exp i-pE) n (E) dE, (7) 

[Q (E) is the state density of the system) reveals two coexisting phases with 
different energies at the same temperature (ferromagnetic and paramagnetic phases). 
The localization of the peaks is determined from those intersection points of the 
microcanonical caloric curve dS/dE where d'^S/dE^ < with the horizontal lines 
representing here the inverse temperature /3 of the thermostat (the ordinary equilibrium 
condition between the thermostat and the system temperature). Since no one peak 
accesses the region where d'^S/dE'^ > at any thermostat temperature, macrostates 
with Ch < are practically inaccessible (unstable) within the canonical ensemble when 
the system of interest is large enough. Such a "forbidden" region is the origin of the 
sudden jump of the canonical average energy (E) at the neighborhood of the critical 
point f3c, which tell us about the existence of a latent heat ql for the conversion of one 
phase into another. 

The close relation between macrostates with negative heat capacities and first- 
order phase transitions in finite short-range interacting systems clarifies that such 
an "anomalous" behavior is far from an unusual feature within the thermostatistical 
description, including also all those systems which have been traditionally considered 
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within the standard Thermo dynamiciJJ]. The anomalous character of such macrostates 
reflects that they are physically admissible within the microcanonical description, while 
they are thermodynamically unstable within the canonical description; indicating thus 
the inequivalence of these statistical ensembles for finite systems. 

As already evidenced, macrostates with C < do not receive a correct treatment 
by the fluctuation relations derived from the usual equilibrium situations of the 
standard Statistical Mechanics and Thermodynamics: in fact, the condition C < in 
Eq.(l2]) cannot be realized in a canonical description, where f3 = constant, or within 
a microcanonical framework, where E = constant so 6E vanishes. Precisely, the 
fundamental aim of this work is to obtain a suitable generalization of thermodynamic 
relation ([2]) in which accounts for appropriately the existence of macrostates with C < 0. 

The generalized fluctuation relation resulting from the present analysis points 
out the role of macrostates with C < in the experimental determination of 
the microcanonical caloric curve (which implies a simultaneous measurement of the 
temperature and energy of a given system), contributing in this way to a re-examination 
of an old question of the Thermo-Statistics Theory: Could it be possible the existence of 
certain kind of complementarity between the energy and the temperaturel Such an idea 
was suggested by Bohr and Heisenberg in the early days of the Quantum Mechanics 
|17] . and so far, it have not received a general consensus in the scientific literature |18j . 

This paper is organized into sections, as follows: First, in section 2, we introduce 
concepts like the generalized thermostat, effective inverse temperature and fluctuation 
in the Gaussian approximation to derive our fundamental result, Eq.( l33i) . Second, in 
section 3, we discuss some implications on themodynamic control and measurements. 
Furthermore, we perform a generalization of this result, Eq. (|52l) . which constitutes 
an "uncertainty relation" in Thermodynamics. Finally, in section 4, we make some 
concluding remarks about several theoretical and practical implications of the present 
formalism. 



2. Extending fluctuation relation 

2.1. Generalized thermostat 

Derivation of a generalized fluctuation relation ([2]) to deal with macrostates with C < 
demands to face the problem of the ensemble inequivalence, to do this, we need to start 
from a more general equilibrium situation of the kind "system-surrounding", where 
macrostates with negative heat capacities can be arisen as thermodynamically stable 
under the external influence imposed by the system surrounding. 

II The consideration of the thermodynamic hmit N ^ oo with N/V fixed in the thermo-statistical 
description of short-range interacting systems is a useful and convenient ideahzation which dismisses 
the occurrence of boundary effects. In practice, any physical system is conforms by a very large but 
finite number of constituents. In fact, the existence of backbending of the microcanonical caloric curve 
as in FIGOis just a finite size effect associated to the presence of interphases during the occurrence of 
the first-order phase transitions (see in [15]). 



Thermodynamic fluctuation relation for temperature and energy 



6 



In thermodynamic equilibrium, the underlying physical conditions of the natural 
environment lead to a situation that a given short-range interacting system obeys the 
Gibbs canonical ensemble ([T]). The natural environment is a good example of a Gibbs 
thermostat, since its heat capacity is so large that it can be considered to be infinite in 
every practical situation. According to Lynden-Bell in ref.[Tl], a system with a negative 
heat capacity Ci < can reach thermodynamic equilibrium with a second system with 
C2 > when < C2 < \Ci\. Obviously, such a condition does not hold when the first 
system is under the influence of the natural environment, since C2 — > 00, and therefore, 
it is not possible to access the macrostates with negative heat capacities. As a corollary 
of this reasoning, the direct observation of the macrostates with negative heat capacities 
of a given system demands that the external influence of the natural environment be 
suppressed. Thus, the system of interest can be isolated (microcanonical ensemble) or 
it can be under the influence of a second system that remains of finite size. The latter 
possibility allows to express the energy distribution function of the first subsystem as 
follows: 



where ^Ic{Et) = Jq ^ ^2 {Et — Ei) ^li (Ei) dEi and as before, ^2 (-^2) is the state 
density of the second system acting as a "finite thermostat", and Ex is the total energy 
of the closed system Et = Ei + E2. 

After reading the analysis presented in the subsequent subsections, it can be realized 
that the consideration of the ansatz ([8]) is sufficient to arrive at a generalized expression 
of the fluctuation relation ([2]): a closed system composed of two independent finite 
subsystems with an additive total energy. However, such a physical picture is only 
admissible when these subsystems are coupled by the incidence of short-range interacting 
forces or when long-range forces are confined to each subsystem which, however, is in 
general nonphysical. Even in this case this assumption presupposes to dismiss the energy 
contribution involved in their mutual interactions Vint- Although this is a licit and useful 
approximation in standard applications of Statistical Mechanics and Thermodynamics, 
it may be unrealistic in the case of mesoscopic systems with short-range interactions, 
and worse, this approach is not applicable in the case of long-range interacting systems. 
The latter case constitutes a typical scenario where macrostates with negative heat 
capacities naturally appear. Obviously, the additivity of the total energy is no longer 
applicable since the interaction energy Vint cannot be considered as a "boundary effect" 
as in the case of large systems with short-range interactions, and often, separability of 
a closed system into several parts is a hypothesis that should be carefully applied [19] . 

Nevertheless, we can find some systems in Nature where it is still possible to assume 
certain separability of a closed system into subsystems despite the presence of long- 
range interactions. Good examples are galaxies and their clusters. Of course, we are 
unable to dismiss the interactional energy Vint in this scenario, but it is reasonable 
to assume as a consequence of the separabihty that this energy contribution could 
be approximated by certain functional dependence of the internal energies of each 




(8) 
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subsystems, Vmt — Vint{Ei, E2)- Thus, the usual additivity of the total energy Ej^ 
could be substituted by the following ansatz: 

= $ (^1, E2) ^E^ + E2 + Vint [El, E2) . (9) 

The number of macrostates of the whole system (Et) is given by: 

fie (Et) = j 6[Et-<^ (^1, ^2)] ^^1 (^1) (^2) dE^dE^, (10) 
which allows to express the energy distribution of the first subsystem as follows: 

p (Er, Et) dE^ = u {E,; Et) fii (^1) dE,, (11) 
and the probabilistic weight uj {Ei, Et) by: 

^ (^1; Et) = 77-^ ! 5[Et-^ [E,, E2)] Q2 {E2) dE2. (12) 

Eg. fill I) constitutes a more general expression than Eq.([H]), providing a better treatment 
for a nonlinear energy interchange between the subsystems as a consequence of the 
nonadditivity of the total energy (AE'i 7^ —/S.E2). This kind of consideration could be 
applicable to both: the case of large systems with long-range interactions and mesoscopic 
systems with short-range interactions. 

Previous discussions have suggested that a general way to account for the energy 
distribution function p{E)dE, which is associated to a general "system-surrounding" 
equilibrium situation, is provided by the ansatz: 

p (E) dE = u (E) n (E) dE, (13) 

where Q (E) represents the state density of the system with internal energy E, and 
uj{E), a generic probabilistic weight characterizing the energetic interchange of this 
system with its surrounding. The above hypothesis is very economical since it demands 
merely: (1) the existence of some kind of separability between the system and its 
surrounding, (2) and those all external influences on the system are fully described 
by the probabilistic weight u (E) . In this work, we are admitting the validity of Eq. (fT3ll 
without mattering the internal structure of the surrounding, and even, the features of 
its internal equilibrium conditions. This last idea is very important and deserves to be 
clarified. 

It is almost a rule that a large system with long-range interactions, initially far from 
the equilibrium, reaches rapidly a metastable equilibrium spending a long time in it [9]. 
If this is the case, the energy interchange of this large system acting as "surrounding" 
of another system cannot be dealt by the expressions ([8]) or ffTTl) . However, this 
metastability does not forbid the applicability of the ansatz ( |T3|) in many physical 
situations. For example, it is well known that the collisionless dynamics of astrophysical 
systems leads to a metastable state where the one-particle distribution function / (g, p) 
depends only on the particle energy e {q,p) = ^p'^ + nup (g), where / (g,p) = f [£ {q,p)], 
whose mathematical form is not Boltzmannian and it is determined from the initial 
conditions of dynamics ^U\. Since the macroscopic behavior of the system in this 
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metastable state is also ruled by its energy, we can expect that this physical quantity also 
rules its interaction with other systems. The admittance of the validity of the ansatz ([T3l) 
for systems in metastable conditions accounts for the claims of some recent authors about 
the existence of non-Boltzmannian energy distribution functions outside the equilibrium, 
overall, in systems with a complex microscopic dynamics, e.g.: astrophysical systems 
[20] and turbulent fluids [21]. A unifying framework for many of these distribution 
functions is provided by the so-called " Superstatistics" , a theory recently proposed by 
C. Beck and E.G.D Cohen [22j, where the presence of a non-Boltzmannian weight u (E) 
seems to be originated from an effective incidence of a fluctuating inverse temperature /3 
at the microscopic level (e.g. on a Brownian particle) obeying the distribution function 
/ (/?) as follows: 

uj{E) = I exp{-pE)fiP)dp. (14) 

Since the microscopic origin and the specific mathematical form of the probabilistic 
weight CO (E) are arbitrary, it is expected that the conclusions, which are derived from 
the ansatz ( fT3l) . may be applicable to a wide range of equilibrium or meta- equilibrium 
situations. In this formula, u (E) constitutes a generic extension of the usual canonical 
weight: 

uj,{E;(3)=Cexp{-[3E), (15) 

which rules the energy interchange of a Gibbs thermostat (a very large short-range 
interacting system). In an analogous way, we shall hereafter refer to the system 
surrounding associated with the weight u (E) as a generalized thermostat. Let us now 
show that this generalized thermostat can be also characterized by certain effective 
inverse temperature /?^, which controls, as usual, the energetic interchange of this 
thermostat with the system under study. 



2.2. Effective inverse temperature 

A straightforward way to arrive at this important thermodynamic quantity is by using 
the probabilistic weight uj (E) in the Metropolis Monte Carlo simulation of the system 
dynamics in this equilibrium condition. The acceptance probability for a Metropolis 
move p{E\ E + AE) is given by: 

p(£|£ + A£) = mm{i;^<§^,l}. (16) 

By assuming that the system size is large enough, consequently, the amount of energy 
|Aii^| <^ \E\, we are able to introduce the approximation: 

and rephrase the acceptance probability (fT6!) as follows: 

p{E\E + AE) ~min{exp[-/3^(E)AE],l}. (18) 
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The evident analogy of this last results with the canonical case allows to define the 
quantity Pui{E) 

[UE) = JIh^ (19) 

as the effective inverse temperature of the generalized thermostat. 

It is possible to verify that this definition is not arbitrary. Besides the obvious 
case of the canonical weight f[T^ where (3i_^ [E) is the inverse temperature of the Gibbs 
thermostat, (E) = P, the present definition drops to the usual microcanonical inverse 
temperature of the system acting as a "finite thermostat" in Eq.([H]), where ut (E) = 
Q2 {E2) {Et) and f3^ (E) = (3^ (^2) = (^2) /dE^, being ^2 (^2) = In [^2 (^2) SEq] 
and E2 = Et- E. 

This inverse temperature is "effective" because it is not always equivalent to the 
ordinary interpretation of this concept in the standard Statistical Mechanics. In order 
to verify this fact, let us consider the equilibrium situation accounted for by Eq. (|TT]) . 
The internal energy of the second subsystem E2 can be expressed as follows: 

E2 = Et-E,-Vm{Ei,E2), (20) 

where the recursive substitution of this last equation in Vmt {Ei, E2) leads to the certain 
nonlinear dependence of E2 on the internal energy of the first subsystem Ei. 

E2 = E2 {El) = Et — (yEi] Ex) ■ (21) 

Thus, the probabilistic weight of Eq.( fT2l) is similar to the case of the "finite thermostat" : 

..(p.p x _ ^2 [^2 (^1)] . . 

u^{Ei,Er)-^^^^, (22) 
but the effective inverse temperature is given by: 

P.iE.) = ^^HE^), (23) 

where the factor z/ (Ei) = —dE2/dEi = dQ {Ei, Et) /dEi 7^ 1 accounts for the existence 
of a nonlinear energy interchange as a consequence of the nonadditivity of the total 
energy. It is remarkable that the "effective inverse temperature" of a very large system 
surrounding depends on the internal energy Ei despite its "microcanonical inverse 
temperature" (32 = 882/ dE2 remains practically unaltered by the underlying energy 
interchange, as in the case of the Gibbs thermostat. 

The physical meaning of the effective inverse temperature (3^ is even more unclear in 
the case where the probabilistic weight u (E) is associated with a system surrounding in 
a metastable equilibrium. Obviously, this kind of inverse temperature has nothing to do 
with the inverse temperature /3 of the integral representation ( fT^ of the Superstatistics 
formalism: while a whole set of values of (3 for each energy E exists, there is only one 
value of (3^ for each value of E. The importance of this new concept relies on the 
possibility to consider a wide class of equilibrium or meta-equilibrium situations in a 
unifying framework where it could be possible to extend the validity of some known 
thermostatistical results. 
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A fundamental identity is the condition of thermal equilibrium. This condition 
commonly follows from the analysis of the most probable macrostate E: 

max{uj{E)n{E)} (24) 

E 

which leads directly to the stationary equation: 

(E) = (E) , (25) 

where (3s (E) is just the microcanonical inverse temperature of the system: 

dSjE) 

Ps [E) = , (26) 

where, as in the introductory section, we take 5* {E) = In W{E). We have demonstrated 
that the quantity jSuj also obeys the ordinary form of the Zeroth Principle of 
Thermodynamics . By analyzing the expression of shown in Eq.f l23l) associated 
with the equilibrium between two separable subsystems with a nonadditive total 
energy, Eq.(l9]), we notice that it is precisely and not (32, which becomes equal 
to microcanonical inverse temperature of the first subsystem (3i during the thermal 
equilibrium, Puj = Pi ^ P2 Pi- This result is essentially the same as the one presented 
by Johal in the ref . |23j . 

Let us now summarize the fundamental properties of the effective inverse 
temperature Puj: 

A. It characterizes the energy interchange "system-surrounding" during the 
equilibrium, in accordance with Eq.( JT8ll . 

B. This concept permits to extend the thermal equilibrium condition, Eq.f l25l) . to a 
wide class of physical situations. 

C. In general, this thermodynamical quantity depends on the internal energy of the 
system E, P^ {E). 

Properties A and B correspond to the ordinary physical notion of temperature. 
A new remarkable property is the general dependence of Pi^ on the internal energy E 
of the system (the property C). This fact indicates clearly that the underlying energy 
interchange not only imposes fluctuations on the internal energy E, but also provokes 
the existence of correlated fluctuations between the internal energy E and the effective 
inverse temperature of the thermostat P^, {SP^SE) ^ 0. This simple property has 
been systematically disregarded by the use of the Gibbs canonical ensemble in standard 
Statistical Mechanics. We shall show in the next subsection that the property C is 
precisely the fundamental key for arriving at a suitable generalization of the fluctuation 
relation ([21). 

Before the end of this subsection, it is important to remark that the probabilistic 
weight: 

uJm{E;Eo) = C6{E-Eo), (27) 
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Figure 3. Typical behavior of the effective inverse temperature is depicted for: a) 
f3t^ fixed for all values of E (canonical ensemble), b) an arbitrary dependence Pai (E) 
(generalized thermostat), and c) E fixed for all values of /3i^ (microcanonical ensemble). 
The black thick point represents the macrostate of thermal equilibrium where Pui — Ps, 
and 0s (E) — dS (E) /dE, the microcanonical inverse temperature of the system. 



which is associated with the energetic isolation or the microcanonical ensemble, can be 
also considered as a limiting case of a generalized thermostat. The application of the 
definition (fT9l) leads to an indeterminate value of /^t^: 



which means that the thermostat inverse temperature f3i^ admits any value when the 
energy of the system is fixed at Eq. This idea is schematically illustrated in FIGJHl 

2.3. Fluctuations in the Gaussian approximation 

We shall suppose in the present analysis that the system and bath are large enough in 
order to support a Gaussian approximation of the energy fluctuations around the most 
probable macrostate E. In addition, we also assume that the energy dependence of 
the effective inverse temperature /J^^ [E) allows for the existence of only one intersection 
point E with the microcanonical caloric curve of the system (3^ [E) in the thermal 
equilibrium condition ( l25l) . This latter requirement is just the condition of ensemble 
equivalence. The average square dispersion of the internal energy {5E'^): 






(29) 



can be estimated as follows: 




(30) 
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The fluctuation 5/3^ = {E)—P^ (i?) of tlie effective inverse temperature can be related 
to the fluctuation 6E of the internal energy of the system: 

dHnuj(E) , ^ 

= dE^^^- (31) 

If we rewrite Eq.f l30l) in the following manner, 

dHnuj(E) . d^S(E) ^ ^ 

and combine Eq. fl3T|) and fl32|) . we arrive at the correlation between the effective inverse 
temperatures of the generalized thermostat and the internal energy of the system as 
follows: 



{5PJE) = 1 + {5E-) . (33) 

This latter identity is a generalized expression of the fluctuation relation ([2]). This fact 
can be noticed by rephrasing fl33|) as follows: 

C = (5^ {5E^) + C {6(3JE) (34) 

by using the microcanonical definition of the heat capacity, Eq.(jl]), being P = Ps (E). 



3. Discussions 



We remark that the fluctuation relation fl33l) accounts for the specific mathematical 
form of the probabilistic weight u (E) in an implicit way throughout the effective inverse 
temperature Once more, this result evidences the fact that the system-surrounding 
energy interchange is effectively controlled by The imposition of the restriction 
SPuj = associated with the Gibbs canonical ensemble ([T]) into Eq. (lMll leads to the 
usual identity C = i^E'^)- However, the restriction 5/3^ = is not compatible 
with the existence of energetic regions with negative heat capacities C < 0. Since the 
microcanonical entropy is locally convex d'^S (-E) /dE^ > in such anomalous regions, 
the identity (!33l) leads here to the inequality: 

(SPJE) > 1. (35) 

This means that any attempt to impose the canonical condition —>■ 0, within regions 
with C < 0, leads to the occurrence of very large energy fluctuations SE oo, 
and conversely, any attempt to impose the microcanonical condition SE — > is 
accompanied by very large fluctuations of the effective inverse temperature of the system 
surrounding SP^j —>■ oo. Remarkably, this behavior suggests the existence of some kind 
of complementarity between the internal energy E and the effective inverse temperature 
Pi^ of the surrounding, which is quite analogous to the complementarity between a 
coordinate q and its conjugated momentum p in Quantum Mechanics! As well, the 
divergence 6Pi^ — » oo when 6E — > is not only applicable to regions where C < 0: any 
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attempt to reduce the energy fluctuations to zero, 6E 0, in Eg. (15^ for any fixed E 
leads to the following result: 

lim (S(3JE) = 1, (36) 

SE^O 

indicating thus the divergence of the effective inverse temperature fluctuations — > oo 
at this limit. The present results allow to conclude: 

(i) While macrostates with C > are accessible within the canonical ensemble, where 
takes place the restriction 5/5^ = 0, the anomalous macrostates where C < can 
be only accessed when 5/5^ 7^ 0, that is, by using a generalized thermostat whose 
energy dependence of its effective inverse temperature /^^^ (E) ensures the validity 
of the inequality (!35l) . 

(ii) The total energy E of the system and the effective inverse temperature of the 
generalized thermostat behave as complementary thermodynamic quantities 
within the regions with C < 0. 

(iii) The imposition of the microcanonical restriction 5E —>■ leads at any internal 
energy E to an indetermination of the effective inverse temperature of the system 
surrounding SjSuj — > 00 as a consequence of the suppression of the underlying 
energetic interchange. 

A deeper understanding of the physical meaning of the previous conclusions is 
reached by discussing their implications on the two standard ways where the external 
influence of the surrounding on the thermodynamic state of the system is used within 
the Thermodynamics: 

• as a control apparatus (thermostat), or 

• as a measure apparatus (thermometer). 

3.1. Implications on the thermodynamic control 

It is well-known that the (inverse) temperature is, in general, a good control parameter 
for the internal energy of a system: the contact of a thermostat with a given constant 
value (3 leads to the existence of small fluctuations of the internal energy 6E oc V^, 
where N is the system size. The remarkable exception is that during the first-order 
phase transitions. Here, a small variation on j3 is able to provoke a sudden change in 
the expectation value of the internal energy (E) of the system due to the multimodal 
character of the energy distribution function in the neighborhood of the critical inverse 
temperature f3c, which provokes very large energy fluctuations 6E oc A^. This physical 
situation was already illustrated in FIGj2]of the introductory section for the case of the 
thermodynamical description of the g = 10 states Potts model on the square lattice 
([6]), where the origin of this anomaly relies on the existence of inaccessible or unstable 
energetic regions with C < 0. 

From the perspective of the thermodynamic control, the fluctuation relation of 
Eq. (l33l) describes the necessary conditions where the external influence imposed by the 




Figure 4. Metropolis Monte Carlo calculations for the 10-state Potts model by using 
Eq. p^ . Consideration of a thermostat with a variable inverse temperature Pui{E) 
can ensure the existence of a unique stable stationary point of the microcanonical 
caloric curve Ps (E) where P^j = Ps, in this way the system can be forced to access its 
anomalous macrostates with C < 0, and at the same time, to eliminate the bimodal 
character of the energy distribution function in the neighborhood of the critical inverse 
temperature /3c- 



surrounding ensures the thermodynamic stability of the system, allowing thus an effective 
control of the internal energy E of the controlled system. The term "effective control" 
means that the internal energy of the system is preserved until the precision of small 
fluctuations around the average value E. This kind of fluctuating behavior is ensured 
by the ensemble equivalence or the existence of only one sharp peak in the energy 
distribution function p (E) dE, which is mathematically expressed by the existence of 
only one intersection point E in the condition of thermal equilibrium (125|) . 

The Conclusion (i) claims basically that macrostates with C < can be forced 
to be thermodynamically stable or accessible by using a generalized thermostat with 
an appropriate energy dependence in its effective inverse temperature Puj (E). This 
situation is clearly illustrated in FIGJH where this particular external influence has 
automatically eliminated the bimodal character of the energy distribution function of 
the g = 10 states Potts model in the neighborhood of the critical temperature shown in 
FIGH 

The divergence SE oo when Sf3uj within the region with C < indicated 
in the Conclusion (ii) accounts for the well-known fact that such macrostates turn 
unstable or inaccessible within the canonical description. This divergence is just a 
consequence of the Gaussian approximation which has been employed to obtain the 
fluctuation relation fl33|) . since the energy fluctuations are actually on the order of 
the system size, SE oc A^, which diverge only in the thermodynamic limit oo. 
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The Conclusion (iii) is just the mathematical result previously obtained in Eg. fl28l) 
and illustrated in FIG|3l Conveniently, the practical implications of this result will be 
discussed in the next subsection. 

As already shown, the consideration of a generalized thermostat with an appropriate 
energy dependence of its effective inverse temperature [E) is a simple but an effective 
consideration to overcome the difficulties associated with the existence of macrostates 
with C < 0. In principle, the practical implementation of the present ideas should 
lead to the development of some experimental techniques which could be particularly 
useful to deal with thermodynamical description of mesoscopic systems which are now 
of interest of Nanosciences, and where the presence of macrostates with C < is not 
unusual. 

Another important framework of applications is to develop more efficient Monte 
Carlo methods to deal with the difficulties associated to the presence of first-order phase 
transitions [21]. A simple and general way to account for the presence of a generalized 
thermostat with a given effective inverse temperature {E) is by using the Metropolis 
method based on the acceptance probability of Eq. ffTSl) . In fact, the microcanonical 
caloric curve (3s i^) of the g = 10 states Potts model shown in FIG|2] and FIGH] were 
obtained using this methodology by calculating the averages (e) and (Puj), where the 
validity of the Gaussian approximation for N large enough and the condition of thermal 
equilibrium (J25l) ensure the applicability of the relations: 

e = {e) and {3^ {e) = {(3^) . (37) 

Numerical results derived from this algorithm agree very well with the ones obtained 
using other Monte Carlo methods [25] . 

The idea of using the generalized thermostat with an appropriate effective inverse 
temperature Puj (E) can be easily extended to other Monte Carlo methods based on 
the Gibbs canonical ensemble. A specific example is the enhancement of the well-known 
Swendsen-Wang (SW) cluster algorithm [16], which suffers from the supercritical slowing 
dowT^ in its application to the q = W states Potts model on the square lattice and is 
unable to capture the C < regime of the microcanonical caloric curve. The thermostat 
inverse temperature (3 enters into this method by the probability for the cluster 
formation p {j3) = 1 — exp (— which is used to generate a new system configuration X. 
While the parameter f3 remains constant in the original SW algorithm, our modification 
consists in the substitution of this parameter by the effective inverse temperature of the 
generalized thermostat which is redefined in each Monte Carlo step, 
The parameter used to generate the configuration takes the value of the 

effective inverse temperature corresponding to the total energy Ei = H]^ (X*) of the 
previous system configuration X\ jS^ = (5^ {Ei). Here, is the Hamiltonian of Eq.(l6]) 
and X = {cTfc} the spin variables. The microcanonical caloric curve is determined by 
using the same relations (|371) of the Metropolis method explained before. In FIGJ5] 

% An exponential divergence of the correlation time tn in the thermodynamic limit iV — > oo, 
Tjv oc exp (TV). 
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Figure 5. The ordinary Swendsen-Wang cluster algoritlim (CE) is unable to account 
for the C < regime of the microcanonical caloric curve /3s (e) = ds (e) /de of 
the q = 10 states Potts model on the square lattice as shown here by our Monte 
Carlo calculations. This difhculty is overcome by the Metropolis algorithm and the 
Swendsen-Wang cluster algorithm based on the use of an effective inverse temperature 
f3uj (GCE), where one can note a remarkable agreement between these algorithms. 
In addition, we also show here the energy per particle e dependence of the second 
derivative of the entropy per particle, k (e) = (9^s (e) /9e^, which has been obtained 
from the fluctuation relation (|33p by calculating the fluctuations quantities (&^) and 
(SfiujSe). This dependence clearly evidence the existence of an anomalous region where 
k(£) > ^ C < 0. 



a comparative study among the ordinary SW method (CE) and the Metropohs and 
SW methods with an appropriate effective inverse temperature (GCE) is depicted. The 
agreement between these last Monte Carlo algorithms and their advantages over the first 
method is remarkable. As clearly evidenced, the present ideas support the development 
of an alternative framework of the well-known multicanonical Monte Carlo methods 
|24j . As in FIGj3]as Monte Carlo calculations shown in FIGHland FIGjSl the following 
effective inverse temperature: 

(3^{E) = r]exp[X{E-Eo)/N] (38) 

was assumed, where t] and A are two real positive parameters controlling the horizontal 
position and the inclination of this dependence respectively. The value A = 
corresponds to the canonical ensemble ([T]), while A —>■ +00 corresponds to the 
microcanonical ensemble (|27|) . More details of these calculations can be found in the 
ref.[26]. 

3.2. Implications on the thermodynamic measurements: An uncertainty relation? 

Bohr and Heisenberg suggested in the past that the thermodynamic quantities of 
temperature and energy are complementary in the same way as the position and the 
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momentum are related in Quantum Mechanics [T7]. Roughly speaking their idea was 
that a definite temperature can be attributed to a system only if it is submerged in a 
heat bath (Gibbs thermostat). Energy fluctuations are unavoidable. On the contrary, a 
definite energy can be only assigned to systems in thermal isolation, thus excluding the 
simultaneous determination of its temperature. Dimensional considerations suggest the 
existence of the following relation: 

6(36E > 1. (39) 

However, these ideas have not reached a general consensus in the literature [ISj.One 
objection is that the mathematical structure of Quantum theories is radically different 
from that of classical physical theories, so that, there are no noncommuting observables 
in Thermodynamics. Interestingly, the result expressed in Eg. (1351) is quite analogous to 
the uncertainty relation fl39|) . 

Before analyzing the implications of the fluctuation relation ( l33l) in the 
thermodynamic measurements, it is important to revise the concept of temperature 
of a system. In the Bohr and Heisenberg arguments explained above, as well as in 
the works of some other authors like the Landau (see last paragraph of pag. 343 in 
ref.[T]) and Mandelbrot [27], a system has a definite temperature when it is put into 
contact with a Gibbs thermostat. In this viewpoint, the concept of the temperature of 
an isolated system is unclear. We think this is a misunderstanding, since the inverse 
temperature f3 that enters in the Gibbs canonical ensemble ([T]) is just the microcanonical 
inverse temperature of the Gibbs thermostat (th), /3 = dSth/dEth- This quantity 
does not properly characterize the thermodynamic state of the system, but instead 
it characterizes the thermodynamic state of the thermostat and its thermal influence 
on the system. Hereafter, we refer to the inverse temperature of a system defined as 
its microcanonical inverse temperature (3s = dS (E) /BE. From this perspective, the 
temperature of an isolated system is a well-defined quantity and does not undergo 
fluctuations, 5E = ^ 5 [3s = 0. 

Although it is possible to obtain the system temperature at a given energy by 
calculating the Boltzmann entropy S (E), the practical determination of the temperature 
is always imprecise. The energy is a quantity with a mechanical significance, and it 
can be determined by performing only one instantaneous measurement on an isolated 
system. The temperature is a quantity with a thermo- statistical significance, whose 
determination demands to appeal to the concept of statistical ensemble. For example, 
it is derived from a statistics of measurements or temporal averages of certain physical 
observables with a mechanical significance usually referred to thermometric quantities 
(e.g.: the length of a mercury column, an electric signal, the pressure of an ideal gas, 
the average form of a particles distribution function, etc.). In practice, the temperature 
(3s of a given system is indirectly measured from its interaction with another system 
(usually smaller than the system under study and referred to thermometer (th)), whose 
internal dependence of its temperature (3th on a thermometric quantity is previously 
known. The fundamental key supporting this procedure is precisely the condition of 
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thermal equilibrium: 

(3s = f3th. (40) 
In our approach, the surrounding can be also used as a measurement apparatus 
(thermometer), so that, the quantity [3th is just the effective inverse temperature [3^. It 
is possible to realize that the generalized thermostat used in the previous Monte Carlo 
calculations has the dual role: to be a control and a measure apparatus at the same 
time. In fact, the energy dependence of the microcanonical inverse temperature j3s of 
the Potts model is a priori unknown, and, we estimate it through the effective inverse 
temperature of the generalized thermostat (3^^ via the condition of thermal equilibrium 
in Eq. (!37l) . This kind of computational procedure exhibits essential features of a real 
determination of the energy-temperature dependence (caloric curve). 

Despite the apparent simplicity, the determination of the energy-temperature 
dependence of a system by using this procedure is rather complex, which is evidenced 
from the two factors affecting its precision: 

A. This measurement process necessarily involves a perturbation on the thermodynamic 
state of the system. 

B. The interaction also affects the thermodynamic state of the thermometer in 
an uncontrollable stochastic way, which provokes the existence of errors in the 
determination of the effective inverse temperature (3^^ used to estimate f3s via the 
condition of thermal equilibrium. 

It is possible to see in the Conclusions (ii) and (iii) that the precision factors A 
and B are rather complementary. The error type A can be characterized in terms of the 
system energy fluctuations 6E, since the inverse temperature fluctuation SPs is directly 
correlated with SE. The error type B is characterized in terms of the effective inverse 
temperature fluctuations 6(3^^, which may also depend on 6E, but in an indirect way. 
According to the Conclusion (iii), any attempt to reduce the perturbation of the system 
to zero, 6E —>■ 0, (the error type A) leads to a progressive increasing of the error type 
B, 6/3^ oo, which affects the estimation of Ps by using (3^. The error type B can 
reduce to zero ^ by imposing the conditions of the Gibbs canonical ensemble (II]), 
which always involves certain perturbations of the system energy 6E ^ (error type 
A). This perturbation is relatively small when the system size is large enough and the 
heat capacity of the system is positive C > 0. However, this situation changes radically 
when the macrostate of the system is characterized by a negative heat capacity C < 0. 
According to the Conclusion (ii), the reduction of the error type B to zero, 5/?^ 0, 
induces the thermodynamic instability of the macrostates with C < and leads to the 
existence of very large energy fluctuations 6E ^ oo. In practice, we should admit the 
simultaneous existence of the errors type A and B, which can reasonably be small and 
unimportant, when the sizes of the system under study and the thermometer are large 
enough. However, such errors are much significant when the system size N is small (in a 
system of few bodies or constituents) that the concept of system temperature becomes 
experimental unobservable, and therefore, physically meaningless. 
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As already evidenced, the fluctuation relation of Eg. (15^ accounts in some way for 
the limit of precision for a determination of the energy-temperature dependence of a given 
system by using a measurement procedure based on thermal equilibrium with another 
system (thermometer). Although this qualitative behavior is quite close to the Bohr 
and Heisenberg intuitive idea of energy-temperature complementarity, the fluctuation 
relation fl33l) does not always support a complementary relationship between the system 
energy E and the effective inverse temperature (3^^ or the system energy E and inverse 
temperature (3s- The reason is that the mathematical structure of Eq.( l33l) does not 
have the form of a complementary relation. Fortunately, this limitation is not difficult 
to overcome. 



3.3. Generalization: The Quantum- Statistical Mechanics analogy 

The derivation of the restricted fundamental result fl33|) relies on the Gaussian 
approximation. However, there is a simple way to overcome this difficulty. The inverse 
temperature fluctuation 6(3s can be expressed for small 6E in the following way: 

6(3s = ^dE; (41) 



therefore, the fluctuation relation ( ]33l) can be rephrased as follows: 

{dvSE) = 1, (42) 

where rj is the difference between the effective inverse temperature of the generalized 
thermostat and the microcanonical inverse temperature of the system: 

V = P.- Ps. (43) 

The validity of result (H2l) does not depend on the Gaussian approximation. In order 
to show this, let us denote by p {E) = uj {E) Q (E) the energy distribution function of 
the ansatz f|T3l) . The distribution function p{E) is not arbitrary. It obeys the following 
general mathematical conditions: 

CI. Existence: The distribution function p{E) is a nonnegative, bounded, and 

different iable function on the set vr C i? of all physically admissible energies E. 
C2. Normalization: This function obeys the normalization condition: 

p (E) dE = 1. (44) 

C3. Boundary conditions: This function vanishes together with its first derivative on 
boundary dn of the set vr: 

WEedn:p{E) = -^p{E)=0. (45) 



The key of the demonstration is the consideration of the following identity: 

_d_ 

"dE 



-^}piE) = r^piE), (46) 
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which allows to associate the quantity rj to an "operator" 17: 

By using the identity (H6|l as well as the properties (C2) and (C3) of the distribution 
function p (E), it is easy to obtain the following results: 

(r^) = TIP (E) dE= I T)p (E) dE = 0, (48) 

J 7T J 77 

{Eri) = / Erip{E)dE= / Ef}p{E)dE = 1. (49) 

J IT J 'K 

The result (1481) is just the condition of thermal equilibrium, Eq. (l25l) . The fundamental 
result (H2|) is immediately obtained from Eq. (H9|) by operating as follows: 

(r^) = ^ I^Et]) - {E) {t]) = {6E6r]) = 1. (50) 

The result fl42p can be rephrased one more time by using the well-known Schwarz 
inequality: 

(M') {6B^) > {6A6Bf (51) 
arriving finally at a definitive result: 

AriAE > 1, (52) 



where Ar] = /(V)" and AE = V7^^. 

Undoubtedly, the expression ( 152|) has the form of "an uncertainty relation" that 
exhibits now a very general validity. It clarifies that the complementary relation actually 
exists for the system energy E and the inverse temperature difference rj between the 
system and its surrounding (acting as a thermostat or a thermometer). Any attempt to 
perform an exact determination of the system temperature f3s throughout the thermal 
equilibrium, At] 0, involves a strong perturbation on the system energy AE 00, 
thus, E becomes indeterminable; any attempt to reduce this perturbation to zero, 
AE —>■ 0, makes impossible to determine the system inverse temperature (3s by using 
the condition of thermal equilibrium since Ar] 00. 

Remarkably, the above thermodynamic complementarity between E and r] is 
quite analogous to the complementarity in Quantum theories. Contrary to what was 
preliminarily objected, this complementarity could be also supported in terms of the 
noncommutativity of mathematical operators E = E and fj = —Be'- 

E,fi = 1 ^ Ar]AE > 1. (53) 

Thus, the relations of Eq. fl53|) can be considered the Statistical Mechanics counterpart 
of the familiar quantum relations: 

[g, p]=i^ AqAp > 1, (54) 

where H = 1. The formal correspondence is: q ^ E and p ^ rj. This Quantum- 
Statistical Mechanics analogy can be also extended to the properties of the distribution 
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functions: p{q) = \ip {q)f and p{E) = uj{E)^l{E), since the quantum distribution 
function p{q) also obeys the properties (CI), (C2) and (C3). Demanding a httle of 
imagination, the behavior of the energy distribution function in the neighborhood of 
the critical point illustrated in FIG|2] possesses a certain analogy with the tunneling of 
the distribution function p (g) (a wave packet) throughout a classical barrier; that is, 
the phase coexistence phenomenon can be interpreted as a "tunneling" of the energy 
distribution function p (E) throughout the canonically inaccessible region. 

Obviously, such a Quantum-Statistical Mechanics analogy should not be 
overestimated, although it is licit to recognize that it is very interesting: (1) Both 
are physical theories with a statistical mathematical apparatus; (2) While Quantum 
Mechanics is a theory hallmarked by the Ondulatory- Corpuscular dualism. Statistical 
Mechanics is also hallmarked by another kind of dualism since it works simultaneously 
with physical quantities with a purely mechanical significance (e.g.: energy) and physical 
quantities with a purely thermo- statistical significance (e.g.: inverse temperature); (3) 
Thermodynamics is the counterpart theory of Classical Mechanics: while Classical 
Mechanics assumes the simultaneous determination of position q and momentum p when 
/i — >^ 0, Thermodynamics assumes the simultaneous determination of the system energy 
E and its temperature in the thermodynamic limit — 0. 

Our approach to uncertainty relations in Thermodynamics constitutes an 
improvement of the works of Rosenfeld and Scholgl [281 EH], which have been also 
based on Fluctuation theory [1]. These authors derived their respective formulations 
from the consideration of a surrounding in the thermodynamics limit, and hence, this 
physical situation actually corresponds to the Gibbs canonical ensemble or in general the 
Boltzmann-Gibbs distribution [18]. Since this work hypothesis presupposes a concrete 
thermodynamic influence, an experimentalist has no free will to change the fluctuational 
behavior of the system in a given thermodynamic state. 

This important limitation is overcome in our work: the fluctuational behavior 
of the system-surrounding equilibrium is modified by considering a different energy 
dependence of the effective inverse temperature (3^1 (E), which presupposes to use a 
different experimental arrangement. In addition, it is not necessary to appeal to 
Statistical Inference theory to arrive at an uncertainty relation, as in the Mandelbrot 
development [30] , which undergoes the same " free will limitation" of the Rosenfeld and 
Scholgl works explained above since it is applicable only to the framework of the Gibbs 
canonical ensemble. Interestingly, no one of the above approaches says anything about 
the relevance of anomalous macrostates with (7 < on the practical determination of 
the system energy-temperature dependence. 

4. Conclusions 

In an attempt to establish an appropriate framework in order to deal with the existence 
of macrostates with negative heat capacities in terms of a fluctuation theory, we have 
introduced the concept of effective inverse temperature P^, which characterizes the 
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thermodynamic influence of tlie surrounding on tlie system under study, and specially, 
allows to extend the condition of thermal equilibrium to a wide class of equilibrium or 
meta-equilibrium situations. Precisely in terms of this quantity we arrive at a suitable 
generalization of the well-known canonical relation between the heat capacity and the 
energy fluctuation, C = i^^^)- The new identity, Eq. (|33ll . defines a criterion capable 
of detecting macrostates with negative heat capacities of the system under study through 
the correlated fluctuations of the effective inverse temperature (3^ of the surrounding 
(generalized thermostat) and the energy E of the system itself. This last constitutes a 
novel kind of mechanics to perform a more effective thermodynamic control on those 
anomalous macrostates hidden behind the phenomenon of ensemble inequivalence, which 
inspires the developments of new experimental techniques of control to deal with the 
thermodynamic description for mesoscopic systems, like those which are now of interest 
in Nanosciences, as well as the introduction of new Monte Carlo methods to overcome 
the difficulties associated with the presence of first-order phase transitions. 

We think that a significant aspect of this paper is to provide new physical arguments 
supporting the existence of a complementary relation, between the system energy and 
temperature, as postulated by Bohr and Heisenberg in the early days of the Quantum 
Mechanics. Both the fluctuation relation (133|) and its generalization (l52l) . which can 
be regarded as an uncertainty relation, impose limitations to the precise determination 
of the energy-temperature dependence (caloric curve) of the system under study by 
using a measurement procedure based on thermal equilibrium with another system 
(thermometer). While this limitation is unimportant in large enough systems, it discards 
the practical utility of some concepts of Thermodynamics in the context of systems with 
a small number of constituents. Surprisingly, these results suggest the existence of a 
remarkable analogy between the Statistical Mechanics and the Quantum Mechanics. 
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